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Abstract 

We investigate squark and gluino pair production at LHC (CMS) 
with subsequent decays into quarks and LSP for the case of nonuni- 
versal gaugino masses. Visibihty of signal by an excess over SM back- 
ground in (n > 2) jets + EJp^^^ events depends rather strongly on the 
relation between LSP, gluino and squark masses and it decreases with 
the increase of LSP mass. For relatively heavy LSP mass closed to 
squark or gluino masses and for {rriq, rrig) > 1.5 TeV signal is too small 
to be observable. For the case when only third generation squarks and 
LSP are relatively light signal is not visible without b-tagging. 
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1 Introduction 

One of the LHC supergoals is the discovery of the supersymmetry. In par- 
ticular, it is very important to investigate a possibility to discover strongly 
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interacting superparticles (squarks and gluino). In ref.[|I| (see, also references 
[@]) the LHC squark and gluino discovery potential has been investigated 
within the minimal SUGRA-MSSM framework where all sparticle masses 
are determined mainly by two parameters: mg (common squark and slepton 
mass at GUT scale) and mi (common gaugino mass at GUT scale). The sig- 
nature used for the search for squarks and gluino at LHC is (n > 0) leptons 
+ (n > 2) jets + E^^^^ events. The conclusion of ref. is that LHC is able 
to detect squarks and gluino with masses up to (2 - 2.5) TeV. 

In this paper we investigate the squark and gluino discovery potential 
of LHC for the case of nonuniversal gaugino masses. Despite the simplicity 
of the SUGRA-MSSM framework it is a very particular model. The mass 
formulae for sparticles in SUGRA-MSSM model are derived under the as- 
sumption that at GUT scale {Mqut ~ 2 ■ 10^^ GeV) soft supersymmetry 
breaking terms are universal. However, in general, we can expect that real 
sparticle masses can differ in a drastic way from sparticle masses pattern of 
SUGRA-MSSM model due to many reasons, see for instance refs. |^, ||, 0]. 
Therefore, it is more appropriate to investigate the LHC SUSY discovery 
potential in a model-independent way []. 

The cross section for the production of strongly interacting superparticles 

PP 99, Q9, m (1) 

depends on gluino and squark masses. Within SUGRA-MSSM model the 
following approximate relations among sparticle masses take place: 

m?~mQ-|-6mi, (2) 

m^o ^ 0.45m 1, (3) 

m^o ^ m^± ^ 2m^o, (4) 

rrig ^ 2.5mi (5) 

The decays of squarks and gluino depend on the relation among squark and 
gluino masses. For rriq > rrig squarks decay mainly into gluino and quarks 

• q^gq 

^The early version of this study has been pubhshed in rcf. pi. 
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and gluino decays mainly into quark- ant iquark pair and gaugino 
•9^ QQxt 

For rriq < rrig gluino decays mainly into squarks and quarks 

whereas squarks decay mainly into quarks and gaugino 

• ? ^ 

• ? ^ dxi 

The lightest chargino Xi has several leptonic decay modes giving a lepton 
and missing energy: 
three-body decay 

• xt ^Xl + l^ + i'. 

two-body decays 

• xf — ' 1%R + 

• xf ^ i^L + 

^x\ + i^ 

• xt^xl + w^. 

Leptonic decays of X2 give two leptons and missing energy: 
three-body decays 

• X^^X? + ^+^, 
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Xf + P + u, 



+ ly 



two-body decay 




As a result of chargino and second neutralino leptonic decays besides 
classical signature 

• (n > 2) jets plus EJp''' 
signatures 

• > 1) leptons plus {n > 2) jets plus E™''^ 

with leptons and jets in final state arise. As mentioned above, these signa- 
tures have been used in ref. for investigation of LHC(CMS) squark and 
gluino discovery potential within SUGRA-MSSM model, in which gaugino 
masses m^o, m^o are determined mainly by a common gaugino mass mi. 



In our study we consider the general case when the relation between 
m^o and rrig is arbitrary. We investigate the detection supersymmetry using 
classical signature (n > 2) jets plus E'™'^''. Signatures with several leptons in 
final state are more model dependent and besides classical signature leads to 
the highest discovery potential within SUGRA-MSSM model. We find that 
LHC squarks and gluino discovery potential depends rather strongly on the 
relation between Xi, 9 and q masses and it decreases with the increase of the 
LSP mass. 

2 Simulation of detector response 

Our simulations are made at the particle level with parametrized detector 
responses based on a detailed detector simulation. To be concrete our esti- 
mates have been made for the CMS (Compact Muon Solenoid) detector. The 
CMS detector simulation program CMSJET is used. The main aspects 
of the CMSJET relevant to our study are the following. 
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• Charged particles are tracked in a 4 T magnetic field. 90 percent recon- 
struction efficiency per charged track with pt > I GeV within \r]\ < 2.5 
is assumed. 

• The geometrical acceptances for /i and e are \7]\ < 2.4 and 2.5, respec- 
tively. The lepton number is smeared according to parametrizations 
obtained from full GEANT simulations. For a 10 GeV lepton the mo- 
mentum resolution Apt/pt is better than one percent over the full rj 
coverage. For a 100 GeV lepton the resolution becomes ~ (1 — 5) ■ 10~^ 
depending on 77. We have assumed a 90 percent triggering plus recon- 
struction efficiency per lepton within the geometrical acceptance of the 
CMS detector. 

• The electromagnetic calorimeter of CMS extends up to |?7| = 2.61. 
There is a pointing crack in the ECAL barrel/endcap transition region 
between \ri\ = 1.478 — 1.566 (6 ECAL crystals). The hadronic calorime- 
ter covers |?7| < 3. The Very Forward calorimeter extends from I?]! < 3 
to |?7| < 5. Noise terms have been simulated with Gaussian distribu- 
tions and zero suppression cuts have been applied. 

• e/7 and hadron shower development are taken into account by parametriza- 
tion of the lateral and longitudinal proffies of showers. The starting 
point of a shower is fluctuated according to an exponential law. 

• For jet reconstruction we have used a slightly modified UAl Jet Finding 
Algorithm, with a cone size of AR = 0.8 and 25 GeV transverse energy 
threshold on jets. 

3 Backgrounds. SUSY kinematics 

All SUSY processes with full particle spectrum, couplings, production cross 
section and decays are generated with ISAJET 7.32, ISASUSY [0. The 
Standard Model backgrounds are also generated with ISAJET 7.32. 

The following SM processes give the main contribution to the background: 
WZ, ZZ, ti, Wtb, Zbb, hi and QCD (2 ^ 2) processes. 

As it has been mentioned above in our paper we consider only classical 
signature (n > 2) jets plus E^^^^ for squarks and gluino detection. We 
considered 3 different kinematical regions: 
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A. rUg S> rriq 

B. rriq » rrig 

C. rriq ~ rrig, rriq > rrig 



We also considered the case when all sparticles are heavy except the 
third generation sfermions and LSP [Q. For the case A squarks pp 



qq production dominates at LHC. The squark decays into quarks and LSP 
Q ~^ QXi lead to the signature 2 jets plus EJ^'^'^^ . For the case B gluino 
pair production pp gg dominates. The gluino decays g qqXi lead to 
{n > 3) jets plus signature. For the case C both squarks and gluino 

are produced pp qq, gg, qg at the similar rate. Their decays give both 2 
jets and (n > 3) jets events. We considered two types of cuts: 

Cuts a. (n > 2) jets with ETjeti > En, ETjet2 > Et2, EJp''' > Etq. 

Cuts b. (n > 3) jets with ETjeti > ^ri, ETjet2 > Et2, Exjetz > Ets, 
^miss y Eto- Cuts a and b are appropriate for the investigation of the 
kinematical points A and B correspondingly, whereas for the point C both 
cuts a and b are useful. We have calculated SM backgrounds for different 
values ETjeti, ETjet2, Exjets, Eto of the cut parameters. Our results are 
presented in Tables 1-2. 

In this paper we considered the case when all squarks have the same mass 
and m^o,m^± > minirrig, rriq). The last requirement leads to the suppression 
of the events with leptons, only classical signature with (n > 2) jets plus 
j^miss essential. The shape of the squark and gluino differential decay 
width depends rather strongly on the relation among squark, gluino and LSP 
masses. We considered different values of squark and gluino masses. We took 
LSP mass equal to m^o = k x min{mg,mq) with k = |, 0.5, 0.75, 0.9. The 
/c = I corresponds approximately to the standard case with universal gaugino 
masses. 



4 Results 

The results of our calculations are presented in Tables 3-6 and in Figures 
1-6. In estimation of the LHC(CMS) gaugino discovery potential we have 

used the significance determined as 5*12 = \/{Ns + N^) — \/{Nb) which is 
appropriate for the estimation of discovery potential in the case of future 



experiments fill]. We also imposed additional requirement that the ratio of 



signal to background events has to be bigger than 0.5 {S/B > 0.5). For the 
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Table 1: Cuts a and the corresponding value of background events for L = 
10 • 



# of cut 


^rniss ^QgYj 


En [GeV] 


Et2 [GeV] 


Nb 


1 


200 


40 


40 


4995783 


2 


200 


100 


100 


3292494 


3 


200 


100 


150 


3097944 


4 


200 


50 


100 


4478452 


5 


400 


200 


200 


180868 


6 


400 


200 


300 


173889 


7 


400 


100 


200 


247991 


8 


600 


300 


300 


8992 


9 


600 


300 


450 


7771 


10 


600 


150 


300 


17662 


11 


800 


400 


400 


1120 


12 


800 


400 


600 


963 


13 


800 


200 


400 


2708 


14 


1000 


500 


500 


229 


15 


1000 


500 


750 


183 


16 


1000 


250 


500 


616 


17 


1200 


600 


600 


38 


18 


1200 


600 


900 


28 


19 


1200 


300 


600 


115 



comparison we also give the values of often used significance M determined 

Ns ™ 
as = — p and 5cr discovery probability 11 for two values of back- 

^f{Ns + Nb) 

ground Nback- Here A^^ = cr^ ■ L is the average number of signal events and 
Nb = ab ■ L is the average number of background events for a given total 
luminosity L. 

As it follows from our results for fixed values of squark and gluino masses 
the visibility of signal decreases with the increase of the LSP mass. This 
fact has trivial explanation. Indeed, in the rest frame of squark or gluino 
the jets spectrum becomes more soft with the increase of LSP mass. Besides 
in parton model pair produced squarks and gluino are produced with total 
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transverse momentum closed to zero. For high LSP masses partial cancella- 
tion of missing transverse momenta from two LSP particles takes place. The 
fact that with rise of LSP mass the E^^^^ spectrum becomes more soft is 
explicitly seen in Figures 1-3. 

Figures 4-6 demonstrate that within well defined cut the number of signal 
events decreases with increase of the LSP masses that complicates supersym- 
metry detection for LSP masses closed to gluino or squark masses. For the 
case when all sparticles are heavy except third squark generation and LSP 
we have found that signal is too small to be observable for all values of third 
generation squark masses and LSP mass. In this cass possible b-tagging will 
be able to suppress the background an make signal observable |. 

5 Conclusion 

In this paper we have presented the results of the calculations for squark 
and gluino pair production at LHC (CMS) with their subsequent decays into 
jets for the case of nonuniversal gaugino masses. We have found that the 
visibility of signal by an excess over SM background in (n > 2) jets + i?™*^ 
events depends rather strongly on the relation between LSP mass Xi ^^'i 
q, g masses. The visibility of the signal for fixed values of squark and gluino 
masses decreases with the increase of the LSP mass. For relatively heavy 
LSP mass closed to gluino or squark masses and for (m^, rriq) > 1.5 TeV 
signal is too small to be observable. 
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Table 2: Cuts b and the corresponding value of background events for L = 



^ of cut 




El [GeV] 


E2 [GeV] 


Es [GeV] 


Nb 


1 


200 


40 


40 


40 


2953667 


2 


200 


100 


125 


150 


957089 


3 


200 


167 


208 


250 


315594 


4 


200 


233 


292 


350 


104932 


5 


200 


300 


375 


450 


79970 


6 


400 


100 


125 


150 


151076 


7 


400 


167 


208 


250 


20392 


8 


400 


233 


292 


350 


9025 


9 


400 


300 


375 


450 


4113 


10 


600 


100 


125 


150 


8774 


11 


600 


167 


208 


250 


4547 


12 


600 


233 


292 


350 


2599 


13 


600 


300 


375 


450 


1701 


14 


800 


100 


125 


150 


1693 


15 


800 


167 


208 


250 


754 


16 


800 


233 


292 


350 


372 


17 


800 


300 


375 


450 


194 


18 


1000 


100 


125 


150 


425 


19 


1000 


167 


208 


250 


234 


20 


1000 


233 


292 


350 


147 


21 


1000 


300 


375 


450 


59 


22 


1200 


100 


125 


150 


99 


23 


1200 


167 


208 


250 


58 


24 


1200 


233 


292 


350 


31 


25 


1200 


300 


375 


450 


22 
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Table 3: The number of events for cut with maximal value of "significances" 
S2 and 5*12 and the corresponding value of the 5a discovery probability for 
L = 10^pb-\ Case A. 



71 ^ 

Mg 












probability 


or discovery 




cutjf^ 


bignal 


S2 


'-'12 


Nback — Nj, 


Nback = 2 • AT;, 


CjreV 


GeV 










(see lab.ij 




2450 


400 


18 


33. 


4.23 


2.52 


0.667 


0.211 


(3050) 


1200 


17 


21. 


2.73 


1.52 


0.055 


0.006 


2050 


350 


19 


176. 


10.32 


6.33 


1.000 


1.000 


(3000) 


1025 
1500 
1950 


19 


63. 


4.72 


2.62 


0.684 
t<0-5 
t<0.5 


0.198 


1550 


260 


16 


1773. 


36.27 


24.06 


1.000 


1.000 


(2000) 


770 

1162 

1395 


13 


1697. 


25.57 


14.33 


1.000 


1.000 


1050 


175 


10 


18793. 


98.43 


58.03 


1.000 


1.000 


(2000) 


525 


10 


10954. 


64.75 


36.26 


1.000 


1.000 




788 


19 


74. 


5.38 


3.02 


0.900 


0.417 




945 


19 


86. 


6.07 


3.45 


0.983 


0.681 


550 


92 


7 


180270. 


275.47 


156.43 


1.000 


1.000 


(2000) 


225 


10 


20210. 


103.85 


61.71 


1.000 


1.000 




412 


10 


9939. 


59.82 


33.24 


1.000 


1.000 




495 


13 


2560. 


35.27 


20.54 


1.000 


1.000 
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Table 4: The number of events for cut with maximal value of "significances" 
5*2 and S12 and the corresponding value of the 5a discovery probability for 
L = IQ^pb-^ . Case B. 



Mg 












probability 


of discovery 




cut^ 


Signal 




'S'12 


Nback = Nt, 


Nback = 2-Nb 


GeV 


GeV 










(see Tab. II) 




2000 


350 


24 


18. 


2.57 


1.43 


0.031 


0.004 


(2950) 


1000 










t<0-5 




1500 


251 


16 


720. 


21.79 


13.75 


1.000 


1.000 


(1950) 


750 

1125 

1350 


15 


456. 


13.11 


7.32 


1.000 
t<0-5 


1.000 


1000 


175 


7 


10970. 


61.94 


34.29 


1.000 


1.000 


(1950) 


500 


11 


3050. 


34.99 


19.73 


1.000 


1.000 




750 


14 


1271. 


23.34 


13.29 


1.000 


1.000 




900 


14 


1100. 


20.81 


11.70 


1.000 


1.000 


500 


94 


2 


1164159. 


530.81 


300.77 


1.000 


1.000 


(1950) 


250 


7 


30710. 


135.85 


83.26 


1.000 


1.000 
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Table 5: The number of events for cut with maximal value of "significances" 
5*2 and 5*12 and the corresponding value of the 5a discovery probability for 
L = 10^pb~^. Case C, cuts a. 



GeV 


GeV 


cut^ 


Signal 


^2 


5*12 


probability of discovery 
(see Tab. I) 


2400 
(2300) 


400 
1150 

1700 


19 


85. 


6.01 


3.41 


0.980 0.661 
t<0.5 

It < '-'-^ 


2100 
(2000) 


350 

1000 

1500 


19 
19 


191. 
63. 


10.92 
4.72 


6.77 
2.62 


1.000 1.000 
0.684 0.198 
t<0.5 
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Table 6: The number of cut with maximal value of "significances" 5*2 and 
5*12 and corresponding value of probability of discovery. Case C, cuts b. 















probability 


of discovery 




cut^ 


Signal 


^2 


'S'l2 


Nback = Nh 


Nback 


GeV 


GeV 










(see Tab. 2) 




2300 


400 


22 


73. 


5.55 


3.16 


0.943 


0.508 


(2400) 


1150 
1700 


24 


18. 


2.57 


1.43 


0.031 

% < 0.5 


0.004 


2000 


350 


16 


283. 


11.06 


6.30 


1.000 


1.000 


(2100) 


1000 
1500 


22 


51. 


4.15 


2.29 


0.446 


0.088 


1500 


251 


14 


2910. 


42.89 


26.70 


1.000 


1.000 


(1550) 


750 
1125 

1350 


14 


1180. 


22.01 


12.45 


1.000 
t<0.5 
t<0.5 


1.000 


1000 


167 


7 


51260. 


191.50 


124.88 


1.000 


1.000 


(1050) 


500 


7 


24240. 


114.74 


68.46 


1.000 


1.000 




750 


10 


5700. 


47.38 


26.64 


1.000 


1.000 




900 


14 


1460. 


26.00 


15.00 


1.000 


1.000 


500 


84 


1 


4330000. 


1604.4 


980.20 


1.000 


1.000 


(550) 


250 


1 


3456000. 


1365.1 


813.11 


1.000 


1.000 




375 


1 


1794000. 


823.3 


460.28 


1.000 


1.000 




450 


6 


108528. 


213.0 


120.83 


1.000 


1.000 
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M(glu;no)= 1 500 GeV, M(squarks)= 1950 GeV 

M(LSP)= 251 GeV 




400 600 800 1000 

Et miss before any cuts 



1200 1400 
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400 600 800 1000 

Et miss before any cuts 
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300 1000 1200 1400 
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Et miss before any cuts 




300 1000 1200 1400 

GeV 

Et miss before any cuts 



Figure 1: The Et miss distribution for different LSP masses (Mg — 
1500 GeV, Mq = 1950 GeV). Case B. 
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M(glu;no)= 1000 GeV, M(squarks)= 1950 GeV 
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Figure 2: The Et miss distribution for different LSP masses (Mg — 
1000 GeV, Mq = 1950 GeV). Case B. 
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c 80 
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M(glu;no)= 2000 GeV, M(squarks)= 1550 GeV 

1(LSP)= 260 GeV 
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Et miss before any cuts 



1400 
GeV 




300 1000 1200 1400 

GeV 

Et miss before any cuts 




400 600 800 1000 

Et miss before any cuts 
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Figure 3: The Et miss distribution for different LSP masses (Mg — 
2000 GeV, Mq = 1550 GeV). Case A. 
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M(gluino)= 2000 GeV, M(squarks)= 1550 GeV 
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Figure 4: The dependence of the signal events on the LSP mass. Case A, 
cut 16a. 
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M(gluino)= 1 500 GeV, M(squarks)= 1950 GeV 
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Figure 5: The dependence of the signal events on the LSP mass. Case B, cut 
16b. 
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M(gluino)= 1 500 GeV, M(squarks)= 1 550 GeV 
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ure 6: The dependence of the signal events on the LSP mass. Case C, cut 
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